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Abstract 

Stargardt Disease (STGD) is the commonest genetic form of juvenile or early adult onset macular degeneration, which is a 
genetically heterogeneous disease. Molecular diagnosis of STGD remains a challenge in a significant proportion of cases. To 
address this, seven patients from five putative STGD families were recruited. We performed capture next generation 
sequencing (CNGS) of the probands and searched for potentially disease-causing genetic variants in previously identified 
retinal or macular dystrophy genes. Seven disease-causing mutations in ABCA4 and two in PR0M1 were identified by CNGS, 
which provides a confident genetic diagnosis in these five families. We also provided a genetic basis to explain the 
differences among putative STGD due to various mutations in different genes. Meanwhile, we show for the first time that 
compound heterozygous mutations in PROM1 gene could cause cone-rod dystrophy. Our findings support the enormous 
potential of CNGS in putative STGD molecular diagnosis. 
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Introduction 

Stargardt disease (STGD) is the most frequent cause of macular 
degeneration in childhood, with a prevalence of approximately 
1:10000 [1]. It is usually diagnosed within the first two decades of 
life and leads to progressive irreversible loss of central vision, 
delayed dark adaptation and a poor final visual outcome. STGD is 
predominantly inherited as an autosomal recessive trait with 
mutations in ABCA4, also known as ABCR, although an autosomal 
dominant form has been also reported [2] . Rare cases of STGD or 
"Stargardt-like" disease phenotypes have been reported with 
mutations in PROM1, PRPH2, VMD2 (also known as BEST1) 
and ELOVL4, which are involved in various physiological 
pathways that are important for macular function [3]. This 
complex arena of genes and clinical features complicates the 
nomenclature in this field [3]; it is unclear how to classify 
individuals with classic Stargardt phenotype. Classic STGD should 
be restricted to only those cases caused by ABCA4 mutations and 
"Stargardt-like" or juvenile macular dystrophy should be used for 
other genetic etiologies. For the purposes of this study, we classify 
our participants with early-onset macular degeneration as 
"putative STGD" cases. 

Stem cell-based therapy shows great promise for the treatment 
of STGD [4] Accurate molecular diagnosis is therefore essential 



for the selection of patients for clinical trials, and is also crucial for 
prenatal STGD diagnosis. However, the genetic diagnosis of 
individuals with putative STGD is an ongoing challenge because 
of the relatively large sizes of some of the genes involved. ABCA4 
and PROM1 are particularly large containing 50 and 26 exons, 
respectively. VMD2, ELOVL4 and PRPH2 have 8, 6, and 2 exons, 
respectively. 

Furthermore, although biallelic mutations in ABCA4 are found 
in most patients with autosomal recessive STGD, there are studies 
which have shown that mutations in the ABCA4 gene are 
responsible for a wide variety of other retinal dystrophy 
phenotypes, such as cone-rod dystrophy (CRD), and retinitis 
pigmentosa (RP) [5,6]. It has also been proposed that individuals 
carrying mutations in ABCA4 may have a higher risk of developing 
age-related macular degeneration (AMD) [1,7]. We therefore 
sought to investigate whether other retinal disease genes besides 
these reported five genes could lead to putative STGD. 

In this study we initially selected known retinal disease genes as 
a gene capture panel and applied a capture next generation 
sequencing (CNGS) approach to identify genetic defects in seven 
putative STGD patients from five independent families. This 
approach was used to test whether additional retinal disease genes 
could lead to putative STGD. 
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Figure 1. Putative Stargardt disease families and genotype. Pedigrees of the putative STGD families in this study. Filled, half-filled and unfilled 
symbols denotes affected, carrier and unaffected status, respectively. Arrow indicates proband. The mutations originated from paternal or maternal 
allele were defined. 
doi:10.1371/journal.pone.0095528.g001 



Materials and Methods 

Patient Recruitment 

This study conformed to the tenets of the Declaration of 
Helsinki and was approved by the Ethics Committee of Eye 
Hospital, W enzhou Medical University. Written informed consent 
was obtained from all participating individuals or their guardians. 
Patients from families (A, B, C, D and E, Figure 1) were recruited. 
Ophthalmic examination was performed for each patient. 
Electroretinography (ERG) and optical coherence tomography 
(OCT) were performed as routine retinal ophthalmic examination. 
A five ml venous blood sample was drawn into an ethylenediamine 
tetraacetic acid (EDTA) sample tube from every subject. Genomic 
DNA was extracted from peripheral blood leukocytes using the 
standard phenol/ chloroform extraction protocols. 

Targeted Exome lllumina Library Preparation 

Genomic DNA was purified and quantified with Nanodrop 
2000 (Thermo Fisher Scientific, DE). The generation of a targeted 
exome lllumina Library was performed according to the 
manufacturer's protocol (MyGenostics, Beijing, China). A final 
library size of 350-450 bp including adapter sequences was 
selected. 



Disease Genes Enrichment, Sequencing Data Generation 
and Bioinformatics Analysis 

A total of 144 disease genes associated with retinal diseases 
including the five known STGD genes (ABCA4, PROM1, 
PRPH2, VMD2 and ELOVL4) were selected by a gene capture 
strategy, using the GenCap custom enrichment kit (MyGenostics), 
as previously described [8] . The enriched libraries were sequenced 
on an IUurnina Solexa HiSeq 2000 sequencer for paired-end reads 
of 100 bp, and analyzed as described previously [8]. Briefly, using 
the Solexa OA the cutadapt (http://code.google.eom/p/ 
cutadapt/), SOAP aligner, BWA and GATK programs to retrieve 
and align to identify SNPs and insertions or deletions (InDels). 
SNPs and InDels were annotated using the exome-assistant 
program (http://122.228.158.106/exomeassistant). Nonsynon- 
ymous variants were evaluated by three algorithms, SIFT 
(http://sift.jcvi.org/), PolyPhen (http: //genetics. bwh.harvard. 
edu/pph2/) and PANTHER (http://www.pantherdb.org/tools/ 
csnpScoreForm.jsp) as described previously, to determine patho- 
genicity [13]. Multiple sequence alignments were performed using 
ESPript3.0 (http:/ /espript.ibcp.fr/ESPript/ cgi-bin/ESPript.cgi). 
Sequencing data were deposited in NIH Short Read Archive 
(SRP036846). 




Figure 2. Color Fundus photographs of the patients. Color fundus photograph of patient 11:1 from family B and patients 11:1 from family E. Both 
of them showed macular atrophy. 
doi:1 0.1 371 /journal.pone.0095528.g002 
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Table 1. Clinical summaries of individuals in this study. 





Family 


Individual ID 


Gender 


Age (Y) 


Age (Y) at onset BCVA at presentation 


Color-vision problems 


A 


111:1 


M 


21 


10 


20/200 


N/A* 




111:2 


M 


13 


11 


20/200 


+* 


B 


11:1 


M 


34 


10 


20/250 


+ 




11:2 


F 


31 


13 


20/400 


+ 


C 


111:1 


M 


23 


10 


20/200 


+ 


D 


11:1 


F 


13 


6 


20/200 


N/A 


E 


11:1 


M 


25 


17 


20/200 


+ 



M.Male; F.Female; BCVA, Best corrected visual acuity; N/A*, No Answer; +*, Affected. 
doi:1 0.1 371/journal.pone.0095528.t001 



Additional Sequencing 

Targeted amplification of ABCA4 and PROM1 sequences was 
performed using PCR (primer sequences and amplification 
conditions in Table SI in File SI). PCR products were sequenced 
on an ABI PRISM 3730 DNA Sequencer. For mutations, 
nucleotide numbering reflects cDNA numbering with +1 corre- 
sponding to the A of the ATG translation initiation codon in the 
reference sequence, according to journal guidelines (www.hgvs. 
org/mutnomen). The initiation codon is codon 1. 



o 
o 

o - 




Physical locations of the captured regions 



Figure 3. Mean coverage of Capture Next Generation Sequencing. A mean coverage of 200 x over the targeted region was achieved; part of 
the sequencing region was showed. 
doi:1 0.1 371 /journal.pone.0095528.g003 



Results 

Clinical Data 

A total of seven patients (two females and five males, Figure 1) 
from five independent families were recruited for this study. 
Clinical summaries, including visual acuity, age of recruitment, 
gender, and relevant ophthalmological findings are described in 
Table 1 , Figure 2 and Figure S 1 in File S 1 . No night-blindness was 
observed in any of the enrolled patients. We noticed the onset of 
disease in patients from family E was later than that of other 
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Paternal mutation Matermal mutation 

A 

c. 4128+1 G>T c. 5218 OT 

GCAGG TACT TGGG T G G T C 

A 1773 V 




C 



c. 1819 G>A c.1819G>A 
G G G C G G G T T G G G C G G G T T 

G 607 R G 607 R 




D 



c. 3430 OG c. 6095 A>G 

ACTGGTCAG GAACATCTT 



C1140W H 2032 R 




E 



c. 1983+1 G>T c. 730 OT 

TTTGTTGAG T G AC C G A C T 

rt R 244 X 




Figure 4. DNA sequence chromatograms. DNA sequence chro- 
matograms of the affected members with putative Stargardt disease. 
The heterozygous peaks of the mutations were pointed out by red 
arrow. The origin of each mutation was being traced from paternal or 
maternal allele. 

doi:1 0.1 371 /journal.pone.0095528.g004 

patients, which led us to suspect that it has a distinct genetic 
etiology. 

Mutation Analysis 

We first calculated the CNGS results for quality. On average, a 
mean coverage of 200 x over the targeted region was achieved 
(Figure 3). Manual checking of sequencing depth of known STGD 



genes (ABCA4, PROM I, PRPH2, VMD2 and ELOVL4) (Table S2 in 
File S 1 ) showed that a mean coverage of 2 1 6 x was obtained. We 
observed that missing coverage of one exon each in PROM1 
(Exon24) and VMD2 (Exon4) genes. We used Sanger sequencing 
for these two missing exons (primer sequences shown in Table S 1 
in File SI). 

Given that autosomal recessive STGD is largely caused by 
homozygous or compound heterozygous mutations, we initially 
scanned the five reported STGD genes for mutations in our 
dataset. Sanger sequencing was used for validation of any 
pathogenic mutations identified, and the segregation of mutations 
was tested in the familial cases. For any missense variants 
identified, computational prediction by three algorithms was used 
to confirm the number of candidate mutations. 

The mutations identified in these cases are summarized in 
Table 2. Briefly, seven mutations (one homozygous) in ABCA4 and 
two mutations in PROM1 were successfully identified in the STGD 
families via CNGS, Sanger sequencing, and co-segregation 
analysis (Figure 4 and Figure 1, Figure S2 in File SI and Table 
SI in File SI). Furthermore, multiple sequence alignments were 
performed and we found that missense mutations in ABCA4 were 
located within a phylogenetically conserved region (Figure 5). 

With the exception of family C all families had compound 
heterozygous mutations. Additionally, because we can define 
whether the mutation originated from paternal or maternal allele, 
we could trace the origin of each mutation (Figure 1, Figure 4). No 
de novo mutations were identified in these families. 

We searched for the mutations identified in multiple databases, 
including the 1000 Genome (1000G, http://www.l000genomes. 
org/), ESP6500 (http://evs.gs.washington.edu/EVS/) and 702 
sample in-house exome database as normal controls. The 
candidate mutations were not present or present at extremely 
low frequency in these databases (Table 2). The molecular genetics 
of these families clearly shows that patients with PROM1 mutations 
have later onset disease, compared with patients with ABCA4 (17- 
years-old vs 6 tol3-years-old). 

We found that in our cohort of putative STGD patients that no 
mutations in additional retinal disease genes besides these reported 
five genes were identified. We did observe additional heterozygous 
DNA variants (Table S3 in File SI). 

Post-ldentified-Mutations for the Clinical Diagnosis 

After the mutations were identified in families, we realized that 
re-evaluation of patients was necessary especially for family E 
because there are no reports that compound heterozygous PROM1 
mutations have been identified for putative STGD. We retrieved 
the clinical data from these families and confirmed that the 
diagnosis of patients from family A-D is STGD. In Family E, we 
found retinal vessels of the patients are moderately attenuated 
(Figure 2), and the macular atrophy was very obvious. The 
electroretinogram (ERG) showed that both the cone and rod 
responses were affected; meanwhile, the cone responses were more 
severely affected than rod responses (Figure 6 A). OCT (Optical 
Coherence Tomography) testing showed the photoreceptor 
segments and the retinal pigment epithelium atrophy seriously 
(Figure 6B). Visual field testing showed central scotomas, while the 
periphery was spared (Figure S3 in File SI). The Arden ratios of 
the electro-oculogram (EOG) were 1.2 and 1.5 (Figure S4 in File 
SI), respectively, compared with that of normal (1.8). Taken 
together, based on the clinical manifestations, the final diagnosis of 
patients from family E was cone-rod dystrophy (CRD). This 
highlights that accurate clinical diagnoses should based on all the 
available clinical data because there are substantially overlapping 
phenotypes between STGD and CRD (Table S4 in File SI). Also, 
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Figure 5. Multiple-sequence alignment of ABCA4 from different species. Multiple-sequence alignment in ABCA4 from different species 
revealed that these mutations were located within a highly conserved region. The mutations were highlighted by red box. Majority of the identified 
mutations (4 of 6) were located within the alpha-helix regions, a-helices are displayed as squiggles; p-strands as arrows; strict fi-turns as TT letters. 
doi:1 0.1 371 /journal.pone.0095528.g005 



it demonstrates that the CNGS method is indeed the best method 
to determine the genetic cause of a heterogeneous disease since it is 
unbiased. It is also the case that CNGS method may be helpful for 
the accurate clinical diagnoses of heterogeneous diseases even if 
the researcher does not have access to clinically well characterized 
patients with different forms of retinal disorders. 

Discussion 

STGD is the most common childhood recessively inherited 
macular dystrophy. The first identified disease gene linked with 
STGD was ABCA4 in 1997 by Allikmets et d, [9], and since then, 
many additional mutations have been identified [1-3,7,10]. Here 
we recruited five families and found that four of them have ABCA4 
mutations, which indicates that it is very informative to screen for 
ABCA4 mutations in STGD families. 

So far, the use of a variety of mutation detection techniques for 
STGD such as SSCP (single-strand conformation polymorphism)/ 
heteroduplex analysis, high resolution melting, microarray, direct 
Sanger sequencing and PCR-Next-Generation Sequencing 
(PNGS), and whole exome sequencing (WES) approaches have 
been reported [2,3,6,10-15]. With the exception of PNGS and 
WES these methods are labor intensive or low throughput 
approaches. Although the PNGS method has the advantage of 
high throughput, it may be a challenge to amplify of all the 
reported gene fragments in one tube. The bioinformatics analysis 
of the results from WGS is still challenging for most laboratories 
and the cost may be prohibitive. In contrast, CNGS allows for the 
comprehensive molecular diagnosis of these heterogeneous genetic 
diseases and has the advantages of speed (Exons of 144 disease 
genes sequenced at one time) and is cost-effective (less than 1/40 
cost of Sanger sequencing); here we demonstrated the usefulness of 
this approach. 

We observed missing coverage of some exons from CNGS 
based molecular diagnosis of STGD, which indicates the method 
still has flaws for applications in clinical genetic diagnosis. It has 
been reported that with deep sequencing, coverage of some 
regions will be missing [16,17]. The main reason for this may be 



the PCR step, which has a bias for amplification of GC-rich or 
repeat fragments under normal PCR conditions. We analysed the 
missing coverage of exons in two genes (PROM1 and VMD2), then 
found there is a repeat A sequence in PROM1 and more than 60% 
GC content in the corresponding exon of VMD2. To fill in the 
missing data, we designed specific primers to amplify these 
fragments. The results from the present study also suggest that 
before searching for the disease-associated mutation, it is necessary 
to check the coverage of the targeted sequence, even though no 
mutations were found in the missing coverage region. 

Definition of a "disease-associated" mutation is a difficult task, 
particularly if no simple functional assays to determine the 
phenotypic effects of specific variations are readily available [1]. 
In general we use the following criteria: if the mutation allele 
frequency is over 5% in the general population as identified by 
bioinformatics analysis of multiple databases, we would treat it as a 
non-pathogenic mutation since STGD prevalence of approxi- 
mately 1:10000. We also checked it whether the mutation was 
reported in the literature or is novel. In this study, we identified 
these nine mutations, including five novels and four previously 
reported, which expands the mutation spectrum of ABCA4 and 
PROM1. 

Analysis of disease allele frequency in specific populations is 
important for clinical genetic diagnosis. There are reports that 
"population-specific" ABCA4 alleles, such as p.G863A/delG863, 
are founder mutations in Northern European patients. We 
searched all these nine mutation in the literature and found five 
novels and four previously reported. Among these four reported 
mutations, p.A1773V in ABCA4 was reported as one of the 
founder mutations (up to 17%) in Latin American population [18]; 
p.R2038W mutation in USA, Estonia and South African 
population; p.R602W mutation in USA, South African population 
[2,3,19]; G607R in the German population [20]. Taken together, 
this study confirmed that these four mutations are pathogenic 
mutations and among these four reported mutations, p.A1773V, 
p.R2038W and p.R602W may have higher allele frequencies since 
they were frequently reported in different populations. We 
observed two mutations (p.R2038W and p.G607R) [1,2], which 



PLOS ONE | www.plosone.org 



5 



April 2014 | Volume 9 | Issue 4 | e95528 



Molecular Diagnosis of Stargardt Disease by CNGS 



Table 2. Identified mutations summary. 





Allele 
frequency 


Family 


Gene 


Identified Mutations (Exon) 


Reported 
or novel 


sir I 


PolyPhen 


DAMTUCD 

KM IN i ncn 


1 000G 


C3rD3UU 


ln-house 


A 


ABCA4 


c.5318C>T;p.A1773V (Exon 38) 


Reported 


D $ 


PD $ 


f 5 


0 


0 


0 






c.4128+1 G>T (Exon 27) 


Novel 


N/A 


N/A 


N/A 


0 


0 


0 


B 


ABCA4 


C.61120T; p.R2038W (Exon 44) 


Reported 


D $ 


PD S 


De $ 


0 


0.00008 


0 






C.1804OT; p. R602W (Exon 13) 


Reported 


D $ 


Benign 


De $ 


0 


0 


0 


C 


ABCA4 


c.1819G>A;p.G607R(Exon 13, Homo*) 


Reported 


D $ 


PD $ 


De $ 


0 


0.000077 


0 


D 


ABCA4 


c.6095A>G; p.H2032R (Exon 44) 


Novel 


D $ 


PD $ 


De $ 


0 


0 


0 






c.3420C>G;p.C1140W (Exon 23) 


Novel 


D $ 


PD* 


De $ 


0 


0 


0 


E 


PR0M1 


C.730OT; p.R244X (Exon 6) 


Novel 


N/A 


N/A 


N/A 


0 


0 


0 






c.1 983+1 OT (Exon 18) 


Novel 


N/A 


N/A 


N/A 


0 


0 


0 



D: Damaging; PD, Possibly damaging; T, Tolerated; DE,Deleterious; N/A, No Answer; 

*Homo, Homozygous mutation; , SIFT (http://sift.jcvi.org/); PolyPhen (http://genetics.bwh.harvard.edu/pph2/). 
PANTHER (http://www.pantherdb.org/tools/csnpScoreForm.jsp); 1000G (http://www.1000genomes.org/); 
ESP6500 (http://evs.gs.washington.edu/EVS/); ln-house, in-house exome database. 
doi:1 0.1 371 /journal.pone.0095528.t002 



have extremely low allele frequency (0.000080 and 0.000077, 
respectively), in databases (Table 2), while no allele frequency data 
(just < 1 / 1 404) is available in the Chinese population due to the 
relatively small in-house sample size. As to real allele frequency of 
the mutations identified in this study, further studies are needed. 
This is consistent with the prevalence of STGD (approximately 
1:10000) since the allele frequency of all the mutations identified in 
this study is not detected in our 702 sample in-house exome 
database. We speculate that the p.G607R mutation may have a 
higher allelic frequency, because the patient from family C has 
homozygous mutation of p.G607R and the parents come from 
different regions. One previous study of STGD in the Chinese 
population, screened part of ABCA4 coding sequence (15 exons) 
and identified two relatively common mutations: T1428M and 
R2040X [21]. To further clarify ABCA4 mutation spectrum in the 
Chinese population, further studies of large sample size are still 
needed. 

So far, more than 200 disease-associated ABCA4 variants have 
been identified (http:/ /www.uniprot.org/uniprot/P78363). We 
manfully mapped these mutations to the ABCA4 protein and 
found the majority of the mutations are located at extracelluar and 
intracellular loops (Figure S5 in File SI). There are four mutation- 



rich regions at the protein level, which suggests that they are in a 
key functional region of ABCA4 (Figure S6 in File SI). For genetic 
diagnosis, it is meaningful to scan these mutation-rich exons. 
Therefore, we manfully mapped all the mutations and found five 
exons (3, 13, 22, 29 and 47) have more mutations per length than 
other exons (Figure S7 in File SI). This indicates that these exons 
may be prioritized for the detection of mutations in ABCA4. 

From our five families study, it is clear that patients with 
PROA11 mutations have a later age of onset, compared the patients 
with ABCA4 mutations. This may be a clinically relevant 
observation. We reviewed the literature and found mutations in 
PROA11 can lead to several diseases (Table 3), including Stargardt 
disease (STGD4) [22,23], retinitis pigmentosa (RP41) [24-28], 
autosomal dominant cone-rod dystrophy (CORD 12, CRD) [16], 
macular dystrophy (MCDR2) [29] and autosomal recessive cone- 
rod dystrophy (CRD) [21]. The major difference between these 
four diseases (STGD, CRD, MCDR2 and RP41) is the results of 
full field ERG and fundus appearance, while night-blindness, the 
ages of onset may be helpful for the clinical diagnosis. CRD is a 
panretinal photoreceptor degenerative disorder with predominant 
loss of cone function that affects the macula early in its course, 
while STGD is a progressive bilateral atrophy of the retinal 



A 

|1 OOjj v i 

I RodL 



Cone 



30Hzl 



B wave 

• o — 



60 80 100 120 140 




Figure 6. Electroretinogram and OCT testing of the patient from family E. A. The electroretinogram (ERG) showed cone responses are more 
severely affected than rod responses (Only B wave has a small peak at 108 ms/56.6uv). This is the key clinical characteristics for diagnosis of CRD. B. 
OCT (Optical Coherence Tomography) testing of the patient from family E showed that photoreceptor segments and the retinal pigment epithelium 
were seriously affected. 
doi:1 0.1 371 /journal.pone.0095528.g006 
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pigment epithelium in the macula, with accumulation of a 
lipofuscin-like substance in the retinal pigment epithelium, and a 
reduced foveal cone ERG. Full field ERG is the key test, 
particularly when patients are asymptomatic and show a normal 
fundus at early stages because full field ERG examination can 
distinguish the effect of degree of cone function and rod function. 
In other words, cone function in CRD would be affected earlier or 
more severely than in STGD as measured by the ERG test. 

Combined the comprehensive clinical examination and genetic 
diagnosis, this is the first report, to our knowledge, to show that 
compound heterozygous mutations in PROM1 could lead to CRD. 
This study also demonstrates that genetic testing can help to 
improve the diagnostic accuracy of heterogeneous disease. 

In the present study, we demonstrated that no additional retinal 
disease genes could cause STGD. This may be due to the limited 
sample size of our study; it also suggests that if additional causative 
genes for STGD exist, that they may be present in a relative small 
fraction of cases. Since there is no more available data from 
additional similar reports, the real fraction still needs to be 
investigated. However, we cannot exclude novel genes, beyond the 
scope of our existing knowledge of retinal diseases. It also indicates 
that, it is very meaningful to scan the mutation in ABCA4 gene 
before screening the all retinal disease genes, since to date mainly 
ABCA4 has been the gene underlying this disorder. 

In summary, we have demonstrated the utility of CNGS 
approach to molecular diagnosis of putative STGD in five 
independent families, with successful identification of disease- 
causing mutations, including five novel mutations. The study also 
provides a genetic basis of the differences among putative STGD 
patients due to different mutations in different genes, which is a 
very significant advance in clinical genetic diagnosis of putative 
STGD. We showed that compound heterozygous mutations in 
PROM1 could cause cone-rod dystrophy for the first time. Our 
findings support the enormous potential of CNGS in putative 
STGD molecular diagnosis. With the progress of next generation 
sequencing technology, higher sequencing quality will be provided 
and its cost will dramatically decrease, which are the key bottom- 
necks for the application of CNGS to clinical genetic testing. Here 
we only showed CNGS results from five families, and more studies 
should be performed before its application as a routine clinical 
genetic testing method. 

Supporting Information 

File SI Contains the following files: Figure SI. Color fundus 
photographs of the patients. Color fundus photograph of probands 
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